Density of states near a vortex core in ferromagnetic superconductors: Application to STM 

measurements 
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We investigate numerically the local density of states (LDOS) in the vicinity of a vortex core in a ferromag- 
netic superconductor. Specifically, we investigate how the LDOS is affected by the relative weight of the spin 
bands in terms of the superconducting pairing, and we also examine the effect of different pairing symmetries for 
the superconducting order parameter. Our findings are directly related to scanning tunneling microscopy mea- 
surements and may thus be highly useful to clarify details of the superconducting pairing in recently discovered 
ferromagnetic superconductors. 
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I. INTRODUCTION 

Recently, UCoGe was added to the distinguished list 
of materials (already featuring UGe2 and URhGe) which 
appear to display coexistence of ferromagnetism and 
superconductivityi^il. While ferromagnetism and conven- 
tional superconductivity may be shown to be antagonistic in 
terms of a bulk coexistent stated, several studies have pointed 
out the possibility of a non-unitary, spin-triplet superconduct- 
ing state coexisting with itinerant ferromagnetism^' ^'^'^''^''" . 
The synthesis of two important phenomena in condensed- 
matter physics, ferromagnetism and superconductivity, is not 
only interesting from the point of view of basic research, 
but has also spawned hope of potential applications in low- 
temperature nanotechnology. 

A number of questions arise concerning the nature of the 
coexistence of ferromagnetic and superconducting order. In 
particular, it is crucial to address i) whether the two long- 
range orders are phase-separated or not, ii) whether the mi- 
croscopic coexistence is spatially homogeneous or not, and 
Hi) what the symmetry of the superconducting order param- 
eter is. Concerning the first question, the answer clearly ap- 
pears to be 'yes', since the onset of superconductivity appears 
inside the ferromagnetic part of the phase diagram^. The sec- 
ond question is, however, still open. Some authors have stud- 
ied spatially uniform coexistence of ferromagnetic and super- 
conducting order^ '* ''"^ ", while others have pointed out the 
intriguing possibility of a spontaneously formed vortex lattice 
state'^ '^ '^, due to the internal field. It has been argued''' that 
a key factor with regard to whether such a spontaneous vortex 
phase appears or not is the magnitude of the internal magne- 
tization M. Finally, although the issue of pairing symmetry 
raised in the third question has not been established conclu- 
sively, the most likely option appears to be a non-unitary, spin 
triplet superconducting state, where the spin of the Cooper 
pair couples to the bulk magnetization through a third order 
term ^ i(dk x d^) ■ M in the Ginzburg-Landau free energy. 
Several studies' ^-'^-'^''^-^" have addressed means by which one 
may identify the pairing symmetry of the superconducting or- 
der parameter in a ferromagnetic superconductor, mainly fo- 
cusing on transport properties. 



Clearly, it would be highly desirable to clarify experimen- 
tal signatures of a possible spontaneous vortex lattice-phase 
realized in a feiTomagnetic superconductor. In this work, 
we present numerical results for the local density of states 
(LDOS) in the vicinity of a vortex-core of a ferromagnetic 
superconductor. Our approach is based on the quasiclassi- 
cal theory of superconductivity, and takes into account sev- 
eral crucial factors such as the depletion of the order param- 
eter near the vortex core in addition to self-consistently ob- 
tained magnetic and superconducting order parameters. Our 
results are directly relevant for scanning tunneling microscopy 
(STM) measurements,^' and may be useful to clarify signa- 
tures of the existence of a spontaneously formed vortex lattice 
and also the pairing symmetry of the superconducting order 
parameter 
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FIG. 1: (Color online) Schematic illustration of the model. 

This paper is organized as follows. In Sec. |lll we establish 
the theoretical framework employed in this work. Namely, 
we use the quasiclassical approximation and solve the Eilen- 
berger equation in the vicinity of the vortex core with appro- 
priate boundary conditions. In Sec. |llll we present our results 
for the spatial- and energy-dependence of the local density of 
states near the vortex core. Specifically, we investigate how 
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the relative weight of the spin bands in terms of the super- 
conducting pairing and different pairing symmetries for the 
superconducting order parameter affect the density of states. 
In Sec. |IV]and|Vl we discuss and summarize the main results 
of the paper We will use boldface notation for 2-vectors, . . 
for 4 X 4 matrices, and ... for 2 x 2 matrices. 



II. THEORETICAL FRAMEWORK 

It is generally believed that the pairing symmetry in ferro- 
magnetic superconductors may be classified as a non-unitary, 
spin-triplet state. Our starting point is the quasiclassical 
Eilenberger equation^^ for such a system, which in the clean 
limit reads (see AppendixlAlfor details) 



iVF • V.g^ + [ep3 + M + A(pf), g^] = 0, 



(1) 



where e is the quasiparticle energy measured from the Fermi 
level, Vi? is the Fermi velocity, and [. . .] is a commutator. The 
exchange field h and the superconducting order parameters 
Act are contained in the terms of M — Miag(r3, T3) in addi- 
tion to 



A(pi.) 



A(pf) 



A(p;^) 

-A*(pf) 



At(pj^) 






Ai(PF) 



(2) 



The matrices pi and are defined in the Appendix. The re- 
tarded part of the Green's function, g^, will have the structure 



_ { g{r,PF,e) /(r,pF,£) 



and must satisfy the normalization condition (g^)^ — 1. Due 
to the internal symmetry relations between the components 
of g^, one may parametrize it very conveniently by means 
of a so-called Ricatti-parametrization-^-^"*. In the absence of 
interband-scattering, the Eilenberger equation decouples into 
two 2x2 equations as follows: 

IVF ■\/ga_+[eT3+(ThTo + Aa_{pF),ga]^0, (3) 
where we have introduced 



Act(Pf) = 



Act(pf) 

-a;(pf) 



(4) 



Note that the gap matrix in Eq. (|4|i is a 2 x 2 matrix in particle- 
hole space, while the gap matrix in Eq. (|2|i is a 2 x 2 matrix 
in spin-space. From Eq. (O, one obtains two decoupled dif- 
ferential equations for and ba' 

ivf • Vact + 2aCT£ - a^A;(pF) - Act(pf) = 0, 
iv;^ • V6ct - 26CTe - 6^ Act(pf) - A;(p;^) = 0. (5) 



Note that the above equations do not have any explicit de- 
pendence on the exchange splitting h. As we shall see later, 
the exchange splitting does however enter implicitly through 
the spin-dependent gaps Act. Note that the magnetic vector 
potential A may be incorporated above simply by a shift in 
the quasiparticle energies: e — > e + evi? • A. In a gauge 
that renders the superconducting gaps to be real, one finds 
that eA eA — V$/2, where $ is the superconducting 
phase associated with the broken U(l) symmetry. There- 
fore, the total Doppler shift in the quasiparticle energies is 
e ^ e — errvvp ■ Vg, where the gauge-invariant superfluid 
velocity is = (V$ — 2eA)/(2TO). Below, we keep the 
distribution of the superconducting phase in the order param- 
eter and consider the case with Ginzburg-Landau parameter 
K ^ 1, for which the magnetic vector potential A may be ne- 
glected. This follows since we are considering only one sin- 
gle vortex, i.e. the zero-field limit, such that only gauge-field 
fluctuations around zero could possibly be relevant. However, 
assuming that the superconductors are strongly type-II with 
K ^ 1, gauge-field fluctuations are suppressed ^^'^^ 




FIG. 2: (Color online) The ratio R between the majority- and 
minority-spin gaps as a function of h/ ^ as obtained from a self- 
consistent, mean-field solution [Eq. ||15H. 

In order to solve the above Ricatti-equations, we follow 
closely the procedure of Ref. 23. Let us consider the term with 
• V in more detail. Assume that we have a cylindrically 
symmetric vortex situated at = = with its axis along 
c. The position vector in this coordinate system then reads 
r = raSL + Tbh. Assuming that the transport of quasiparticles 
primarily takes place in the a — b-plane, we may define the 
Fermi velocity as 



vf — fF(cos Ok + sin 9h) = ufv 



(6) 



and its orthogonal vector u = — sin Ok + cos Oh. Thus, the 
position vector r may also be expressed as r = xv + yu, 
where we have defined 

X — ra COS 9 + Tb sin 0, y = —Va sin + ri, cos 0. (7) 

Using the new coordinate system v — u, the Ricatti equations 
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may be rewritten as 

ivpdxaa + [2e 
wpd^b^ - [2e 



= 0, 
0, 



(8) 



where Ua = ac,{x,y) and Ag. — Acr{x,y). The above 
equations may be solved by imposing boundary conditions 
for {aa,ba} in the bulk of the superconductor. The Ricatti- 
equations with e > for and b^y are stable for integration 
from X (— oo) and a; — > oo, respectively (opposite for 
£ < 0);^ The boundary conditions then read: 



acr[x (-oo)] 

ba{x — > oo) 



|A„ 



(e-v/e2-|A.P)/A, 



(9) 



The superconducting order parameter A^ is now modelled in 
the presence of a vortex centered at ra — rt, = 0. In general, 
the superconducting order parameter may be written as2i 



A,(r,0,£) = A,^oX^(0,£)i^(r)e'™^ 



(10) 



assuming a vorticity m. Here, Aq is the gap magnitude, 
x{d, s) is a symmetry factor for the gap (taking into account 
both anisotropicity and frequency-dependence), F{r) models 
the spatial depletion of the gap near the vortex core, while 
tan (j) — ri,/ra- We will here restrict our attention to an even- 
frequency, p-wave symmetry, which is believed to be the most 
likely candidate for the order parameter in ferromagnetic su- 
perconductors. Assuming that the angular symmetry is the 
same for both the majority and minority spin gaps and consid- 
ering the usual case of m = 1, we explicitly have 



A„{r,e) = A^^oxl^*) tanh ( 



\/x^ + y^\ X + ly 



(11) 



In what follows, we will compare the cases x(^) = cos 9 and 
x{S) — e'^, and also investigate how the LDOS changes de- 
pending the relative weight of the superconducting instability 
in both spin-bands. The normalized LDOS for spin species a 
is given by 

N,{r,e)^ — Re{(l-a,6,)/(l + a,6,)}, (12) 
Jo 27r 

and we introduce the total LDOS in the standard way as 



iV(r,e) = ;^iV,(r,e)/2. 



(13) 



To account for a finite quasiparticle lifetime r, we let e 
e + i5 where S ~ t~^. From now on, we fix S = 0.lA|.o 
and comment further upon the role of inelastic scattering in 
Appendix Ib] 

Even if the exchange field h is absent from the Eilenberger 
equation, the LDOS is not independent of the value of h. The 
reason for this is that the magnitude of the superconducting 
gaps depend on the strength of the exchange splitting. Follow- 
ing the approach of Refs. 8. 9, we derive from a weak-coupling 



mean-field theory that the self-consistent solution of bulk su- 
perconducting gaps in the T ^ limit may be written as 



Acr,o/t^o = cexp[-l/(.gv/l + crh/^i)], 



(14) 



where the prefactor is equal to c ~ 2.43 for a pa;-wave sym- 
metry [x{Q) = cosO], c = 2.00 for a chiral p-wave symme- 



ti-y [x{e) 



Here, 
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VqNq is the weak-coupling 



constant which we set to g ^ 0.2 and loq is the typical fre- 
quency width around Fermi level for the bosons responsible 
for the superconducting pairing. Above, Vq is the strength of 
the pairing interaction, iVo is the LDOS at Fermi level in the 
normal-state and /i denotes the Fermi energy. The reader may 
consult Appendix ICl for a derivation of Eq. ( fT4b . We find that 
the ratio between the majority- and minority-spin gaps may be 
written as 



A 



T,0 _ 



A 



= R{h/ fi) = cxp 



(15) 



when assuming that h/fi E [0, 1) (shown in Fig. |2]l. In UGe2, 
the energy splitting between the majority and minority spin 
bands was estimated^ to lie around 70 meV, which yields R ~ 
1.42 when assuming fJ. — I eV. 



j?j--wa\ii hyiiiiuGtry fv(fi) = cos 5] 




FIG. 3: (Color online) Normalized LDOS in the vortex core for a 
Pa;-wave symmetry [x(^) ~ cos6] using several values of h/ /i. 



III. RESULTS 

We begin by plotting the energy-resolved LDOS in the vor- 
tex core (va — ri, — 0) for an order parameter which has 
line nodes in momentum space. Such an order parameter was 
recently proposed to be realized in UGe2 by Harada et al.^, 
and it was moreover argued that the superconducting pairing 
only took place in the majority spin-band. To investigate how 
the relative magnitudes of the majority and minority spin gaps 
affect the LDOS in the vortex-core, we plot the LDOS for sev- 
eral values of the ratio h/ ji in Fig. [3] As usual, the LDOS is 
strongly enhanced for subgap values due to the existence of 
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FIG. 4: (Color online) Normalized LDOS in the vortex core for a 
chiral p-wave symmetry [x(^) ~ e'^] using several values of h/fi. 




FIG. 5: (Color online) Total and spin-resolved LDOS in the vortex 
core for the a) p^^-wave symmetry symmetry [x(^) = cos6] and 
the b) chiral p-wave symmetry [x{9) = e'*]. We have here used 
h/fi = 0.15. 



bound states within the vortex cora^. The presence of two 
gaps in the system should manifest itself in the form of non- 
monotonous behaviour in the subgap spectrum, but it is not 
possible to discern such behaviour unambiguously from Fig. 
[3] This effect may be masked by strong inelastic scattering, 
modelled here by the parameter 5, which effectively smears 



the LDOS. The effect of increasing the exchange field is seen 
to suppress the deviation from the normal-state LDOS. This 
may be understood by noting that the minority-spin gap is 
strongly reduced with increasing exchange field, and that the 
corresponding increase of the majority-spin gap is not able to 
compensate for the suppressed regime of bound states within 
the core. 

We next study the chiral p-wave symmetry analogous to the 
A2 -phase in liquid "^He, and plot the energy-resolved LDOS 
for several values of h/ fi in Fig. |4] Although the qualitative 
behaviour is quite similar to Fig. [3] there are two important 
distinctions. First, one notices that the chiral symmetry ap- 
pears to have a much more pronounced influence on the LDOS 
quantitatively, yielding a larger zero energy-peak and larger 
subgap dips. This is in fact opposite what one would have ex- 
pected from tunneling conductance measurements of pj.-wave 
and chiral p-wave superconductors, respectively. For such 
measurements, the zero-energy peak becomes much larger in 
the pj.-wave case than in the chiral p-wave case. Secondly, the 
subgap features associated with the presence of two gaps are 
enhanced in Fig. |4]compared to Fig. [3] The non-monotonous 
behaviour for subgap energies is present for all curves in Fig. 
m but the features indicative of multiple gaps are most clearly 
seen for h/fi — 0.15, manifested through an additional inflec- 
tion point before the normal state LDOS is recovered. These 
differences could be helpful in discriminating between differ- 
ent types of pairing symmetries in ferromagnetic supercon- 
ductors. 

In order to show more clearly the contribution from each 
spin-band to the LDOS near the vortex core, consider Fig. |5] 
where we plot the total LDOS and the contribution from each 
spin band for a) x(6') — cos 6 and b) x(6') = e'^. The rise 
of the LDOS following the gap edge Acr.o of each spin band 
occurs at different energies due to the exchange splitting. This 
is revealed in the total LDOS as kinks located at two distinct 
energies, which offers the opportunity to obtain explicit in- 
formation about the relative magnitude of the two gaps. The 
qualitative features are the same in Fig. |5]a) and b), but they 
are quantitatively more pronounced in the chiral p-wave sym- 
metry case. This may be due to the fact that the chiral p-wave 
gap has a constant magnitude (|x(6') | = 1), while the pj;-wave 
gap varies in magnitude upon traversing around the Fermi sur- 
face. Therefore, the LDOS is more strongly affected in the 
chiral p-wave case. 

We now study the resolution of the LDOS in real space for 
a fixed energy in Fig. |6l We have chosen R = 2, correspond- 
ing to h/fi ~ 0.14 and also chosen the line node symmetry 
x{d) — cos 6*. In all cases, the plots in Fig. |6]display a two- 
fold spatial symmetry, in accordance with the superconduct- 
ing order parameter. The zero-energy peak present for 
e = evolves into a dip-structure at the vortex core upon 
increasing the quasiparticle energy. The deviation from the 
normal-state LDOS is still significant even at distances ^ 2^ 
away from the vortex core around e/A| o = 0.5. The qual- 
itative features are the same for the chiral p-wave symmetry 
in Fig. |7] although the symmetry is now circular due to the 
isotropy of the magnitude of the gap (|x(^)l = 1) 
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FIG. 6: (Color online) Normalized LDOS in the vicinity of the vortex core at three different quasiparticle energies, using R — 2 with a 
Px-wa.\e symmetry [x(^) ~ cos 6]. A two-fold symmetry is observed in agreement with the symmetry of the order parameter. 




IV. DISCUSSION 



In our calculations, we have chosen a real gauge for both su- 
perconducting order parameters A^r, cr =T, i- If the two spin- 
bands are completely independent, there is no phase-locking 
between the order parameters which fixes the relative phase 
Av = i/| — vi, where v^r is phase associated with the bro- 
ken U(l) symmetry. The existence of two such phases would 
imply that a U(l)xU(l) symmetry is broken in a ferromag- 
netic superconductor, and would in principle allow for two 
critical temperatures which may differ in magnitude. How- 
ever, if the two spin-bands do communicate by means of e.g. 
spin-orbit coupling or impurity scattering, a term of the form 
— A cos(A!^) will appear in the free energy describing the sys- 



tem. This corresponds to a phase-locking scenario where the 
sign of A determines whether Ai/ = or Ai/ = tt is the 
energetically preferred relative phase. Above, we have decou- 
pled the two spin-bands such that the relative phase of A| and 
A^ is of no consequence. Taking into account scattering be- 
tween the spin-bands would require solving coupled Ricatti- 
equations and investigating the effect of phase-locking explic- 
itly, which is beyond the scope of this paper. 

In Fig. |5j we plotted the relative contribution from the two 
spin-bands to the LDOS near the vortex core to clarify how the 
LDOS may give decisive clues about whether both spin bands 
partake in the superconducting pairing or not. In principle, it 
might be possible to probe explicitly the spin-resolved LDOS 
by using a strong ferromagnetic STM and contrasting parallel 
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and antiparallel relative configuration of the exchange fields in 
the FMSC and the ferromagnetic STM tip. The experimental 
realization of this particular proposal is nevertheless probably 
challenging. 

V. SUMMARY 

In summary, we have numerically studied the local density 
of states (LDOS) in the vicinity of a vortex core in a ferro- 
magnetic superconductor. Specifically, we have investigated 
what influence the exchange field and the symmetry of the su- 
perconducting order parameter exhibit on both the spatially 
resolved and energy-resolved LDOS. The symmetry of the 
spatially resolved LDOS near the vortex core as revealed by 
STM-measurements should give decisive clues about the or- 
bital symmetry of the superconducting order parameter^2i22i22 
while the energy-resolved LDOS could provide important in- 
formation about the presence of multiple gaps in the system. 
Our results should be comparable to experimentally obtained 
data, both qualitatively and quantitatively, and may thus be 
helpful in clarifying the nature of the superconducting pairing 
in ferromagnetic superconductors. 
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APPENDIX A: MATRICES AND QUASICLASSICAL 
THEORY 



for nice reviews. We follow here closely the notation of 
Ref.— . Our starting point is the following Hamiltonian: 

H = J2 [ dr^l(r,i)(- ^1-/173) >(r,t) 

^ — ' J \ ZTTl / ap 

drdr'[A.(r,r')^t(r)^t(r) 

+ AUr,r')Mr')Mr)]- (A2) 

The Heisenberg equation of motion for the above Hamiltonian 
is obtained in the standard way: 

i9t/53*(r,i) - J dr'7?(r,r',i)*(r,t), 
H{r, r', t) = i{r)S{r - r') - A(r, r'), ^V) = 

Zm 

A(r,r') =diag[AT(r,r'),Ax(r,r')]. (A3) 

For simplicity, we consider only the retarded component of 
the Green's function in what follows, since the system is 
specified exclusively by in an equilibrium situation. It is 
defined as 

G^^(l,2) = -Mh - t2)([^a(l), V^;^(2)]+>, (A4) 

where the notation (1, 2) refers to the spatial and time coor- 
dinates: (1) = (ri;fi). We explicitly write the '+' sign as 
a subscript to denote an anticommutator; it is else implicitly 
understood that the notation [...] denotes a usual commutator. 
Similarly, the anomalous Green's function is given by 



The matrices used in this paper are defined ai 



,31 



F^^(l,2) = -ie(ti-<2)([^a(l),^M2)]+>. 



(A5) 



Pi 



n 

Tl 



T2 



P2 



-1 

1 



1 

1 

1) ' ^ " (0 1 



1 

-1 



T,; 



-iri 

ITl 



n 



P3 



1 Q 

-1 



(Al) 



Let us briefly sketch the way to obtian the quasiclassical 
Eilenberger equations for a non-unitary, spin-triplet super- 
conducting state coexisting with ferromagnetism. For further 
background information on the quasiclassical theory of super- 
conductivity, the reader may consult e.g. Refs. l32H33ll34ll35ll36l 



One may construct 4x4 matrices in combined particle-hole 
and spin space, known as Nambu space, in the following man- 
ner: 



Note that G(l,2) is a generalized Gor'kov Green's func- 
tion, which contains information about processes occuring at 
length scales comparable to the Fermi wavelength. Such in- 
formation is lost upon applying the quasiclassical approxima- 
tion. Using the Heisenberg equation of motion Eq. iA3h we 
obtain 
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= %<5(l-2). 



(A7) 



To arrive at Eq. ([T]i, it is convenient to introduce the mixed 
representation which shifts the frame of reference to a center- 
of-mass system. We define 

R = (ri + r2)/2, r = ri - r2, 
T={h+t2)/2, t = ti-t2, (A8) 

such that 



G'^(l,2) = G^{K+-,T 



t r t 



The Fourier-transformation of Eq. ( |A9l l yields 

G'*(p,R;T,£) = / dre-'P'- / dte'*'^G"*(l, 2). 



(A9) 



(AlO) 



An exact solution for G^(p, R; T, e) is very hard to achieve, 
but the situation is considerably simplified if one is willing to 
neglect all atomic-scale fine structure effects that are included 
in G^. These give rise to a rapidly oscillating part in the so- 
lution for G^, and rewriting the Green's function through Eq. 
( lAlOb allows us to integrate out this unnecessary information 
(at least for our purposes). This approximation may be ex- 
pected to yield satisfactory results if the energy of the phys- 
ical quantities involved in the problem, e.g. exchange field 
and superconducting order parameter, are much smaller than 
the Fermi energy. Assuming that only particles in the vicinity 
of Fermi level will take part in physical processes, one only 
needs to retain the direction of the momentum at Fermi level 
in the p coordinate. 

As this Appendix is only meant as background informa- 
tion for the Eilenberger equation, we do not show all the 
details leading from Eq. (IA7b to Eq. ([TJ here. The cal- 
culations are nevertheless fairly straight-forward, and con- 
sist of first switching to a mixed representation, then Fourier- 
transforming the variables, and finally performing the quasi- 
classical approximation 

/dCpG^ ^p = |^. (All) 
TT J 2m 



APPENDIX B: INELASTIC SCATTERING 

The choice of 5 = 0.lA|^o is motivated by the fact that the 
zero-energy peaks observed in experiments are usually limited 
from above to roughly a factor of five times the normal-state 
value of the LDOS, which we reproduce with this particular 
choice of 5. Choosing 5 smaller (corresponding to a longer 
quasiparticle lifetime since 5 = r^^) causes the zero-energy 
peak to grow substantially, as shown in Fig. [8] In general, the 
inelastic scattering rate does not have to be proportional to the 
gap at all and our choice of 6 = O.lAj o is simply chosen to 
compare the scattering rate against a familiar quantity. 
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FIG. 8: (Color online) Normalized LDOS in the vicinity of the vortex 
core at three different values of the inelastic scattering rate 5 = r~^, 
using apa;-wave symmetry with h/ /i = 0.15. 



APPENDIX C: DERIVATION OF EQ. 

The gap equation may be obtained by starting out with a 
Hamiltonian assuming a non-unitary triplet pairing state co- 
existing with itinerant ferromagnetism^''^, namely 



INAP 



k 



kcrcr^kcrcr 



kcr 



Ckcr Akcrcr 
^kcrcr ~Ck(T 



Cko- 

'-kc 



(CI) 



Here, / is the ferromagnetic exchange coupling constant, N is 
the number of lattice sites, M denotes the magnetic order pa- 
rameter (dimensionless), while fekcra is the Cooper pair expec- 
tation value. Diagonalization of this Hamiltonian produces: 



kcr 



INAP 



(C2) 



kcr 



where {7kCT,7k(j} are new fermion operators and the eigen- 
values read 



Eka - \ S,] 



'ko" 



(C3) 



Above, ^k is the kinetic energy measured from Fermi level. 
By minimizing the free energy, one obtains the gap equation 
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for the superconducting order parameterji 

Akaa = -^y] ^kk'..||^tanh(/3^k'./2). (C4) 

k' 

Assuming that the gap is fixed on the Fermi surface in the 
weak-coupling hmit, one may write in general 

v„As,o') = -VoY-{e)[Y''{e')]*- (C5) 

where (9) are basis functions for the angular dependence 
of the interaction. To model p^-wawe and chiral p-wave pair- 
ing, respectively, we use Y°'{9) — — cre"^^ and Y"'{9) = 
cos 9. Conversion to integral gap equations is accomplished 
by means of the identity 

^E/(^k.)= / dsN-ie), (C6) 

k •' 

where N'^{e) is the spin-resolved density of states. In three 
spatial dimensions, this may be calculated from the dispersion 
relation by using the formula 

With the dispersion relation ^kcr = £k ^ cr/M — /i, one obtains 



In their integral form, the gap equation reads 



(C9) 



Consider now T = 0, where the integral may be done analyt- 
ically to yield; 



A,,o = ca;oe-i/3^^+^, a i (CIO) 



where we have defined M = IM/fi = h/ fi, i.e. the exchange 
energy scaled on the Fermi energy. Moreover, c is a numeri- 
cal prefactor which depends on which symmetry one consid- 
ers (pa;-wave or chiral p-wave) while g is the weak-coupling 
constant. The important influence of the magnetization is that 
it modifies the density of states, which affects the supercon- 
ductivity gaps. For M = 1, i.e. an exchange splitting equal to 
the Fermi energy, the minority spin gap is completely sup- 
pressed. Thus, the presence of magnetization reduces the 
available phase space for the minority spin Cooper pairs, sup- 
pressing the gap and the critical temperature compared to the 
pure Bardeen-Cooper-Schrieffer case. 
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